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Abstract 

We discuss transport equations resulting from relativistic diffusions 
in the proper time. We show that a solution of the transport equation 
can be obtained from the solution of the diffusion equation by means of 
an integration over the proper time. We study the stochastic processes 
solving the relativistic diffusion equation and the relativistic transport 
equation. We show that the relativistic transport equation for massive 
particles in the light-cone coordinates and for the massless particles in 
spatial momentum coordinates are related to the (generalized) Bessel dif- 
fusion which has an analytic solution. The solution describes a particle 
moving in a fixed direction whose frequency distribution is the Bessel pro- 
cess. An approach to an equilibrium in a moving frame is discussed. We 
formulate the equilibrating diffusion and transport processes in a Lorentz 
covariant way. 

1 Introduction 

There are many versions of the relativistic diffusion and of the relativistic trans- 
port (see [1][2][3] for reviews). However, if we assume that the diffusing particle 
does not change its mass, the diffusion evolves in the proper time and that it 
is generated by the second order Lorentz invariant differential operator, then 
its definition is unique as shown by Schay [4] and Dudley [5]. The theory can 
be generalized to describe a particle diffusing on a manifold (pseudoriemannian 
background metric) [6] [7]. We have extended these models to a diffusion of a 
particle with spin [8] . We have shown [9] that in the models of the relativistic 
diffusion the energy and the angular momentum grow exponentially fast in the 
proper time. Such a growth is rather unphysical in a steady world around us. In 
refs.[10][9] we have studied a modification of the model by an addition of some 
drag terms (friction) resulting from the requirement that a definite equilibrium 
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state is achieved at a large time. In [11] we have shown that a relativistic dif- 
fusion equilibrating to the Bose-Einstein equilibrium is a linearized version of 
the Kompaneets equation [12] applied for a long time to a description of pho- 
tons diffusing in an electron gas [13]. There are phenomena in astrophysics 
[13] [14] [15] as well as in particle physics [16] [17] which can be described by 
relativistic kinetic and diffusion equations. The general theory of relativistic 
kinetic equations has been developed in [18] [19]. Various forms of relativistic 
kinetic equations can be derived from relativistic quantum field theory [20] [21] 
[22]. These kinetic equations are applied to the quark-gluon plasma [17] and to 
the heavy ion collisions [16] often in the diffusion approximation. 

In our earlier papers [10] [8] we have concentrated on the diffusion in the 
proper time. The proper time is a reasonable mathematical tool (see [23] for its 
use in relativistic dynamics) to keep explicit Lorentz invariance. However, we 
should be able to describe the relativistic phenomena also in the coordinate time 
(we call it the laboratory time). In relativistic dynamics there is an equation 
relating the proper time with the coordinate time and we can exchange one into 
the other. The probability density as a function of the proper time trajectories 
should be independent of the proper time. This requirement is equivalent to an 
evolution equation of the probability density in the coordinate time (the kinetic 
equation). It is natural to interpret an analogous requirement (the diffusion 
probability distribution defined on the Minkowski space-time is independent 
of the proper time) as the kinetic equation (the transport equation) in the 
stochastic case. 

We can obtain the transport equation by a random time change as well. 
This is a stochastic analogue of the change of the proper time (associated with 
a moving particle) into a laboratory time. In this paper the random time change 
(which is a standard notion in the theory of diffusion processes [24]) is formulated 
as a calculational tool useful for analytic as well as numerical calculations. The 
method involves an integration over the proper time analogous to Feynman's 
proper time formulation of quantum mechanics [25]. 

In the case of particles of moderate velocity (non-zero mass) the meaning of 
the proper time diffusion can be elaborated by means of an expansion of mo- 
menta in powers of (mc) _1 . In this paper we study in detail the ultrarelativistic 
case of massless particles (then the proper time is a formal affine parameter on 
the trajectory) and concentrated beams of massive particles when only light- 
cone momenta p + are relevant. We perform probabilistic calculations of the 
integral over the proper time by an application of some results of Yor [26] [27] . 
It comes out that the solution of the transport equation obtained as an integral 
over the proper time from the solution of the relativistic diffusion is expressed 
by Bessel functions. This is not accidental as the transport equation coincides 
with the Bessel diffusion. We compare our probabilistic proper time method of 
the solution of the transport equation with the theory of Bessel diffusions. 

In general, the relativistic diffusion will have no limit for a large proper time 
(no equilibrium) . For the same reason the solution of the transport equation will 
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have no limit as the laboratory time tends to infinity. The diffusion equation 
must be supplemented by some drag terms describing a friction if the diffusion is 
to achieve an equilibrium. The equilibrium distribution depends on the velocity 
of the frame where it is observed. So, the equilibrium is a covariant but not an 
invariant notion. In this paper we formulate the dependence of the relativistic 
dynamics on the equilibration reference frame in a covariant way. 

The plan of the paper is the following. In sec. 2 we review the mathematical 
scheme of the relativistic diffusion theory in the proper time. In this section we 
also formulate and prove the basic theorem which allows to obtain a solution 
of the transport equation if the diffusion in the proper time is known. In sec. 3 
we discuss various coordinates which will be applied in subsequent sections. 
In sec. 4 transport equations for relativistic diffusions with drifts leading to an 
equilibrium are discussed. In sec. 5 the relativistic diffusion process as a function 
of the Brownian motion is obtained in an explicit form. In sec. 6 a diffusion 
of massless particles is discussed. Sec. 7 is the main part of this paper. We 
study a relation between the relativistic diffusion in the proper time and the 
transport equation in the ultrarelativistic case of massless particles and the 
light-cone variables for massive particles. We solve the transport equation by an 
application of the theorem of sec. 2. These equations are solved also by another 
method. The result is that the diffusing particle has the momentum whose 
probability distribution is the Bessel process. In sec. 8 we discuss the frame 
dependence of the equilibrating diffusion process. We formulate the theory in 
terms of relativistic invariant variables. The results are summarized in the last 
section of this paper. 

2 The general scheme of the relativistic diffu- 
sion and relativistic transport 

Following [4] [5] [10] we consider an evolution in the proper time of points on the 
phase space (x,p) 

preserving the mass-shell H+ 

rf v VvV„=^V»=vl-\? =m 2 c 2 . (3) 

In eq.(2) F can be a random force leading to a diffusive behaviour of a test 
particle described by coordinates (x,p). From eq.(l) and (3) (for m > 0) it 
follows that r really has the meaning of the proper time. Then, the evolution 
of a scalar function of the trajectory is determined by the equation 

d T ^ = g^> = (p^d x u + A)cp, (4) 
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where (in the case of a diffusion) A is an 0(3, 1) invariant second order differen- 
tial operator; differentiation over space-time coordinates has an index x whereas 
differentiation without an index concerns momenta. The probability density $ 
evolves according to an adjoint equation 

d T $ T = g*$ T = (- P »d* + A*)*r (5) 

resulting from (x € R 4 , p € R 3 ; fourvectors have Greek indices and threevectors 
Latin indices) 

J dxdp<fi T (x,p)$(x, p) = J dxdp<f)(x, p)<j> T (x, p). (6) 

A is the second order differential operator such that Al = and the quadratic 
term is negatively definite. Then, there exists a diffusion process (x T ,p r ) start- 
ing from (x,p) at r = such that 

4> T = E[cj)(x T ,p T )] =/C T 0, (7) 

where 

/C r (x,p;x',p') = £[£(x r (x,p) - x')6(p T (x,p) - p')]- (8) 

From cq.(6) 

$ T - /c;$ (9) 

where 

K*{x, p; x', p') = E[6(x T (x',p') - x)6(p T (x' , p') - p)]. (10) 

The probability distribution <f> is independent of r if it satisfies the transport 
equation 

G*^ = 0. (11) 

We define an equilibrium distribution &e as an x° independent solution of 
the transport equation. Let 

$ = (12) 
Then, * satisfies a diffusion equation 

T tf = (13) 

where 

5 = -^+^ (14) 

and A = Q^AQe is a second order differential operator closely related to A. 
The solution of eq.(13) can be expressed by a stochastic process 

* T = .E[*(x T ,p T )]. (15) 
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Then, if <I> is a solution of the transport equation (11) then 'J is a solution 
of a diffusion equation ( we set t = ^- > 0) 

«9 t 1> = -pV x * + — AV. (16) 
Po Pa 

A solution of the diffusion equation (16) can again be expressed by a diffusion 
process (x t ,p t ) 

* t = pt)]. (17) 

We can establish a relation between the solutions of the diffusion equations (13) 
and (16) 

Theorem 

Let (x T ,p T ) be the diffusion process (15) solving the proper time diffusion 
equation (13) then the solution of the transport equation (16) with the initial 
condition *(x, p) reads 

*t(x,p) = — / drEh{y: T ,i> T )po(T)d(t- — [ dsp (s))} (18) 
mcj L \ mcJo /J 

where 



Po(s) = \J m 2 c 2 + p^. 

Proof: Let us first check the initial condition. For small t the proper time 
t is also small and /J" po{s)ds ~ poT, then x T ~ x and p T ~ p. An integration 
over r gives ^(x, p) when t — > 0. Next, for a derivation of the diffusion equation 
(16) we apply the formula 5(/(r)) = |/'| _1 <5(t — r(t)) in order to express the 
time t in the diffusion process (x r ,p T ) by t. Then, = Po(t) -1 cancels 

the po term in eq.(18). There remains a function of the process (x T ,p T ) at 
the time r(t) in eq.(18).This is the random time change well-known from the 
theory of diffusion processes. It is proved in [24] that after the random change 
of time eq.(16) is satisfied (we shall still discuss the random time change in 
sec.7,eqs.(76)-(79)). 

We shall apply the formula (18) for an explicit calculation of expectation 
values. Although in general the formula (18) may be difficult to use for an 
analytic treatment it still can be very useful for numerical calculations. 

We can express the solution of the transport equation (16) in terms of a 
kernel JC t 

*t=^t*. (19) 
Then, from eqs.(17)-(18) the kernel can be expressed in the form 

= % V £drE[5{^ p) - xO<J(p T (p) - p')) P o(r)S(t ± J Q T dsp ( S ))} 
= ^(x t (x,p)-x')5(p t (p)-p / )]- 

(20) 
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In this section we have applied a notation suggesting the choice of spatial 
momenta p as coordinates on the mass-shell (3). We shall apply some other co- 
ordinates. The definition of the adjoint in eq.(6) is with respect to the Lebesgue 
measure which is not relativistic invariant. A density factor ensuring the invari- 
ance is contained in <I>. In the next sections we shall apply light-cone coordinates. 
Then, the evolution in xq is replaced by an evolution in x_ = x° — x 3 but the 
integration in eq.(6) is still with respect to a Lebesgue measure. 



3 Various coordinate systems 

The relativistic diffusion is generated by the Laplace-Beltrami operator on the 
mass-shell p 2 = m 2 c 2 . If we choose the spatial momenta p as coordinates on 
the mass-shell then the diffusion generator reads 

2 7 " 2 .4 = A H = (5 jk + m - 2 c- 2 p 3 p k )£-£- + Sm^ 2 ^. (21) 

We shall also use n 2 = mT 2 cr 2 ^ 2 as the diffusion constant. 

The spatial momenta are convenient for a physical interpretation. However, 
for a derivation of explicit solutions some other coordinates are more useful. 
Let us consider the Poincare coordinates (91, 92, 93) on Ti+ which are related to 
momenta p as follows 

Ps+Po = — , (22) 

93 

P3-po = (<?i +q 2 +g 3 )> 

93 

Pi = , (23) 



Pi 



93 

mcq2 



93 

where (73 > 0. Then, the metric is 

ds 2 = (mc) 2 q^ 2 (dql + dq\ + dqf) (24) 

and 

(mc) 2 A H = q 2 (d 2 + 3\ + d 2 ) - q 3 d 3 . (25) 

The Poincare coordinates are closely related to the light-cone coordinates 
(p+,Pa) (where a = 1,2) 

P+=Po+P3- (26) 
We have q a = PaP^ 1 and q 3 — mcp^ 1 . Then 

A H = dl + d\ + {mc)- 2 p\d\ + {mc)- 2 p\d\ 

+ (mc)- 2 p 2 dl + 2{mc)- 2 p+p a d a d + + 2{mc)- 2 p 1 p 2 d 1 d 2 (27) 
+ {mc)- 2 'Sp + d+ + (mc)- 2 '3p a d a 
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where d a = jjp ■ The generator Q in the light-cone coordinates is 

g = p-dl+p + d*- Pa dZ + ^A H . (28) 

The time evolution is in x_ instead of x$. In the formulae of sec. 2 we replace po 
byp + and x by The integration measure in eq.(6) is dx\dxidxj r dx-dp\dpidp+. 

4 Transport equations with friction 

In order to achieve an equilibrium at a large time we add a friction term K to 
the diffusion generator (21) 

K = Ki JL. (29) 
The transport equation (11) with friction reads 

— 2 dp, dp k \°3 k ^ " l C PjPk)^ 2 dp k Pk dpk k 

Eq.(30) determines the drift K if <f>£ is space-time independent 

K fe = KVA^Irfe+m-^-Vjft) - |(mc)- 2 p fe )*£- (31) 

Let us consider a class of solutions <&e depending solely on p . We write in 
the form 

$ E =Po 1 exp(/(/?cp )). (32) 

The Pq 1 factor ensures the Lorentz invariance of the measure dpp$ . From 
eq.(31) we obtain 

- — Pk Pcpof((3cpo). (33) 

We choose $ in the form 

$=p - 1 exp(/( C /3p ))* (34) 
Then, the transport equation (16) reads 

= — (pV x + i)* (35) 
Po 

where 

A = K j 8 i +A (36) 

and .4 is defined in eq.(21). 
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Explicitly 

A* = -PiTT-* 

+^(^ + m- 2 c- 2 PjPk )4-^ + + Ic/J/'^o^oKafr*- 

(37) 

(the formulae (35)-(37) are not obvious but come out from calculations). The 
transport equation (without friction) in light-cone coordinates is 

p_^+P + ^*-p o 5„*-^AJ r * = (38) 

where x + = x° + x 3 and x_ = x° — x 3 . The adjoint of the operator (27) in 
eq.(38) is with respect to the Lebesgue measure dp + dp\dp 2 . 

We look for an equilibrium solution of the transport equation in the form (it 
is not normalizable with respect to the integration over the transverse momenta 

Pa) 

$ £ =p; 1 exp(/(/3cp + )). (39) 

We shall still discuss equilibrium distributions of the form (39) in sec. 8. Here, we 
only point out that the light-cone coordinates are appropriate for particle beams 
moving in the direction of the third axis such that the remaining momenta p_ 
and p a are small. An analogue of eq.(31) gives 

K = K_d+ + K a d a = ^l3cK 2 p\.f'd + + n 2 p + d+ - \^p a d a (40) 

for the friction leading to the equilibrium (39). A particular (non-normalizable 
but relativistic invariant) distribution corresponds to / = 0. We note that 
the change of coordinates (pi,p 2 ,P3) — » (pi,p 2 ,P+) transforms the relativistic 
invariant measure dpidp 2 dp 3 pQ X into 

dp 1 dp 2 d P+ p+ 1 . 

With the friction (40) the diffusion generator (4) reads 

2j- 2 A = d 2 + d\ + {mc)- 2 p\d 2 + + (mc)- 2 p 2 d 2 

+ (mc)- 2 p 2 d 2 + 2{mc)- 2 p + p a d a d+ + 2{mc)- 2 p l p 2 d 1 d 2 (41) 

+5{mc)- 2 p + d+ + f3c T ^p 2 + fd+ 

and the diffusion operator A of eqs.(13)-(14) and eq.(36) in the evolution equa- 
tion for the proper time probability distribution 

T tf = (-p-dl - p + d x _y + Pa d a t> + A)V 

is 

2 7 - 2 i = dl + 81 + {mc)- 2 p\d 2 + + (mc)- 2 p 2 d 2 

+ {mc)- 2 p 2 d 2 + 2(mc)- 2 p + p a dad + + 2{mc)- 2 p 1 p 2 d 1 d 2 (42) 
+ {mc)- 2 p+d + + (mc)- 2 2p a d a + f3c^p 2 + fd+ + 20c^p +Pa fd a . 
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5 Stochastic equations 

In the Poincare coordinates the diffusion process is a solution of the linear 
stochastic differential equations (-f 2 — m 2 c 2 n 2 ) 

dq a = Kq 3 db a , 

a = l,2, 

K 2 

dq 3 = — —q 3 dr + nq 3 db 3 = -K 2 q 3 dr + nq 3 o db 3 , (43) 

where Stratonovitch differentials are denoted by a circle and the Ito stochastic 
differentials without the circle (the notation is the same as in [29]). The Brow- 
nian motion appearing on the rhs of eqs.(43) is defined as the Gaussian process 
with the covariance 

E[b a (T)b c (s)} = 5 ac min(T, s). (44) 

The solution of eq.(43) is 

q 3 (r) = cxp(-K 2 r + Kb 3 (r))q 3 (45) 

and 



(46) 



Qa(r) = q a + K / q 3 (s)db a (s). 
Ja 

In the light-cone coordinates 

P± = Po ± P3 (47) 

we have 

3k 2 

dp + = K 2 p + dr + Kp + o db + = —^-p + dr + np + db + , (48) 
dp a = K 2 p a dr - Kp a db + + jdb a (49) 

where a = 1, 2. 

The solution of eqs.(49) is 

Pa(r) = exp(« 2 r - Kb + (r))p a +7/ exp(«; 2 s - nb + (s))db a (s) = e a + S a . 

_ (50) 

The similarity between the solutions in light-cone coordinates and the Poincare 
coordinates is not accidental. It follows from the relation (23) between these 
coordinates. 

A relativistic transport equation can be important in applications to high- 
energy plasma. In an electromagnetic field there is an additional drift term in 
the diffusion equation (4) 

K=—F jv p v d\ (51) 
mc 
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where denotes the electromagnetic field tensor. The stochastic equations 
are linear and can be solved if the only components of F are F12 = B and 
^30 = E. In such a case the stochastic equations (48)- (49) read 

dpi = n 2 pidT + aBp 2 dr — np\db + + jdbi, (52) 

dp 2 = n 2 p 2 dr — aBpidr — np 2 db + + jdb 2} (53) 
dp + = K 2 p + dr + aEp+dr + np + o db + , (54) 

where 

e 

a = — . 

mc 

Eqs.(52)-(54) should be supplemented by equations determining the coordinates 
dx a — m~ 1 p a dT, dx+ = mT x p-dT and 

dx- = m~ p+dr. 

A particular solution $e of the transport equation in an electric field (51) 
(F30 = E, the remaining F^ v = in eq.(51)) 

-p + d x _<$> - p_<9* $ + p a d x $ + A*$ - aEd+p+Q = 0, 

where A* is the adjoint of A (eq.(41) with / = 0), is 

*s=P+ 1 . (55) 
Let $ = p^ 1 ^ then the proper time diffusion equation for * takes the form 

a T * = (-p+d x _ - p-d% + Pa d x +A + aE P+ d+)y. 

The transport equation reads 

d x ^ = + P+Va^tf + P+ 1 A^ + aEd+H , (56) 

where A is defined in eq.(42) with / = 0. In the light-cone coordinates the 
stochastic equation for the process p+ of eq.(15) in an electric field E (but 
without friction) is 

dp + = ^n 2 p + dr + aEp + dr + np + db + . (57) 

It has the solution 

P+{t) = cxp(aEr + k6+(t)). (58) 

If we have a friction leading to the Juttner equilibrium distribution (without an 
electromagnetic field and /' = —1) then the stochastic equation for p + is 

dp + = ^n 2 p + dr — ^l3cn 2 p 2 + dT + np + db + . (59) 
We discuss its solutions in sec. 7. 
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6 Massless particles 

There is a substantial simplification of stochastic equations if m = 0. In order to 
obtain the limit of zero mass we let Sj k — > in eq.(21). The diffusion generator 
in the limit m — > reads 

A H=P lPk w;w: + 3 P k w:- ( 6 °) 

The transport equation for a diffusion with a friction is (when m=0we use 
only the diffusion constant n 2 ) 

If $e depends only on momenta then from eq.(61) 

= ^(ll^PiPk IPk)*E. (62) 

Then, eq.(13) reads (in the massless case the proper time is just an affine time 
parameter without any physical meaning) 

T tf T = \n 2 p jPk dW k ^ T + 2K 2 Pj di* r + \K 2 p jPk {di ln$ E )d k t> T - p M fi£* T . 

(63) 

The transport equation follows from eqs.(62)-(63). In an explicit form 
— * — —n —^i 4- —n nu — — ^ 

dx" * - p p 3 dx 3 * + 2 Po Pjyk dp . dpk * / g4 -. 

We discuss in more detail the Juttner equilibrium distribution [34] 

®e = Pq 1 cxp(-c/3p )- (65) 
Then, from eq.(62) we obtain 

CK 2 

K k = -—flpkpn. (66) 

The stochastic equation for the diffusion (15) in the proper time with the friction 
(66) reads 

3k 2 cfin 2 
dpj = —^-pjdr —pj\p\dr + tvpjdb. (67) 

From eq.(67) we obtain the stochastic equation for p = |p| 

dpo = ^-Podr - ^-f3pldT + np db. (68) 
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7 Solution of the transport equation 



In this section we wish to apply the Theorem of sec. 2 in order to derive a 
solution of the transport equation. Let us begin with the spatial momenta as 
coordinates on the mass-shell H + . The stochastic equation for the transport 
diffusion (x t , p t ) (17) of a massive particle solving eq.(37) in p coordinates reads 

dpj = m^±c 2 (3f'dt + ^^pjdt + mcs/^p^ndbj 

+v / mcK(p§ — mcp 2 )p~ 2 pjpdb (69) 

= 2UpR c 2f3f dt + ^^p 3 dt + sfcnp-^e^dbn, 

dx j = -cp 1 pQ 1 dt, (70) 
here we defined the square root of the metric (the second order terms in eq.(21)) 

e]e% = S jk + {mc)- 2 PjPk . 

The process (69)-(70) gives the solution of the transport equation (37) in the 
form 

*t(x,p)=£?[*(x t (x,p),p t (p))]. 

Unfortunately, the non- linear equations (69)- (70) are difficult to solve explicitly. 
We could also consider the proper time equations in these coordinates discussed 
extensively in [10]. Then, we could apply the Theorem of sec. 2. However, 
such a method would be fruitful only in numerical calculations. Computing the 
expectation values (18) can be more efficient then solving eqs.(69)-(70) directly. 

The light-cone coordinates are more useful for analytic solutions. The formula 
(18) for the solution of the transport equation in these coordinates takes the 
form 

y(X-;X + ,X a ,p+,Pa) 

= -kJ(T dTE V(Xa(T),X + (T),p+(T),p a (T))6(x- - ± /J" dsp+(s)ds)p+(T) , 

(71) 

where p is the stochastic process (15) generated by A of eq.(42). If there is no 
friction then the stochastic equations of the diffusion (42) are linear and can 
be solved explicitly (see eqs.(50) and (58)). Then, the expectation value of a 
function h of p a has a representation in terms of its Fourier transform 

E[h(pi(T),p 2 (T))} = J dk 1 dk 2 h(k 1 ,k 2 )E[exp(ikaPa(T))] 

= J dkidk 2 h(ki,k2)E exp(i/c a e a (r)) cxp ^ - \k a k a J exp(27 2 s + 2jb + (s))ds^j 

(e a is defined in cq.(50)).The joint distribution law of 6, j exp b, J exp 2b can be 
explicitly calculated [30]. Hence, a solution of the transport equation (including 
the electromagnetic field (51)) can be obtained in a form of an integral kernel. 
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We consider here a simplified version of the Theorem of sec. 2 when the initial 
condition <3/ depends solely on p + . Let us note that p+ enters the generator of the 
diffusion with the factor (mc) _1 (or with k = (mc)^ 1 ^). Hence, it is irrelevant 
for moderate velocities. Our restriction to probabilities depending only on p + 
in fact applies to the ultrarelativistic limit when the transverse momenta can be 
neglected. In such a case the transport equation with the Juttner friction (39) 
(/' = —1) but without the electric field reads 

<9_* = 4p+<9+* + 4< 9 +* - ±c[3K 2 P+ d + t>. (72) 

The transport equation with the electric field E and without friction (/ = in 
eq.(39)) defined in eq.(56) can be written explicitly as 

= + (if + aE)d+^> . (73) 

Eq.(71) gives the solution of eqs. (72)-(73) in the form 

= m~ 1 f™dTE tf(p + (T,p + ))*(s- - £ /J" d*p + (s,p+))p+(T,p+) . 

We can solve the diffusion equations (72)-(73) by means of the stochastic process 
(17). For eq.(72) the stochastic equation is 

K 2 1 

dp + = ~^-dx- — -(3cK 2 p + dx- + K\fp\db (75) 

whereas the stochastic process solving eq.(73) satisfies the equation 

,k 2 



dp + = ( — + aE)dx_ + ny/p + db. (76) 

We shall discuss solutions of eqs. (72) and (73) together with the transport equa- 
tion for massless particles because all these equations are related to the Bessel 
diffusion. 

Before solving the equations let us explain the text-book method [24] of the 
random change of time in eq.(67) (then we show that the Theorem of sec. 2 is 
its efficient realization). Let us treat the formula (1) for x°(r) as a definition of 
the proper time t 



/' 

Jo 



- x da (77) 



(here p s is the solution of the proper time diffusion (67) ; x° (r) can be defined 
implicitly by (77)). We can sec from eq.(77) that t depends only on events 
earlier than x°. As a consequence Pj{x°) = pj(T(x )) is again a Markov process. 
Then, differentiating the momenta and coordinates according to the rules of the 
Ito calculus [29] we obtain the following Langevin equations (for mathematical 
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details of a random change of time sec [24] [29] ; a random time change from x° 
to r is discussed in [28]) 

d Pj (x°) = ^\p\-i pj dx° - ^f Pj dx° + K Pj \p\-*db, (78) 

dx j =-p j \p(x°)\- 1 dx°. (79) 

Let ^(x, p) be an arbitrary function of x and p and (x(x° , x, p), p(x°, p)) the 
solution of eqs.(78)-(79) with the initial condition (x, p) then 

* t (x,p) = E[*(x(a; ,x,p),p(z ,p))] (80) 
is the solution of the transport equation 

^ = + &V*lfelfe* " + (81) 



d 

dx" 



with the initial condition \&(x, p). We can prove that eq.(81) is satisfied by 
differentiation of eq.(80) and application of the rules of the Ito calculus [29] (or 
the well-known relation between Langevin equation and the diffusion equation). 

We return to the diffusion in the proper time in order to perform the average 
over the random time t(x ) as prescribed by eq.(18) (this formula is treating the 
random change of time as a computational tool). From eqs.(67)-(68) it follows 
that 

Pt I Pr | ~ 1 = n = (cos (j> sin 9, sin <\> sin 9, cos 9) = const (82) 

is time independent. Hence, according to eq.(18) the solution of the transport 
equation (64) takes the form 

tt t (x,|p|,n) = — / dT£[tf(x(T,x),|p T |,n)j(t / |p s |ds)|p T |l, (83) 

mcj I V mcj a ) J 

where 

x(r,x)=x-/ ds|p s |n. (84) 
Jo 

In order to calculate the expectation value (83) we apply the results of Yor 
[26] [27] who obtained (using the Feynman-Kac formula for an exponential po- 
tential) the joint distribution of 

B^(T) = b T + HT (85) 

and 

(r) = £ ds exp (2B M (s)) . (86) 

We need some rescalings (using ab(s) = b(a 2 s)) in order to bring our formulae 
to the Yor's form. We write 

cx P (k 2 s + nb(s)) =cxp^2^2k 2 j+6(k 2 ^^ . (87) 
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It follows that n = 2 in Yor's formula. 

After such an rescaling an expectation value of a function of B T = k 2 t+k6(t) 
and A(t) = ds exp 2B S is [26] [27] 



E 



f(b{t),A{t)) ] = exp(-^) / F(v, 4, \p\y) exp(2«) 



(88) 



exp(- -^(l + exp(2w)))% 1 exp(w), ^-)y 1 dydv 
where (v G R,y G i?+) 

0(r,r) = r(2^r)-^cxp(^) / °° exp(-|^ - r cosh(O) sinhtf) sin(^)d£. (89) 

The T-integral in the Theorem can be calculated by an application of another 
Yor's formula 

/•OO 

/ dr C M-^r)0(r,r)=I^(r), (90) 
Jo 

where I a is the modified Bessel function of order a [33]. Hence, 



J °° (ItE[F(B(t),A(t))] = £fF(v,£\p\y)e X p(2v) 
CX P (~ + cx P( 2w ))) / 2(j/ _1 exp(w))y- 1 dydz; 



(91) 



We can apply the formula (91) only to the expectation values (18) of diffusions 
without friction (because only in this case the integrand (18) depends solely 
on A and B) corresponding to the Jiittner distribution with / = (this is 
the Lorentz invariant but non-normalizable distribution dxdp<&E = dxdpp^ 1 ). 
Then, with the proper rescalings we have (we omit the index fi = 2 in A and B) 



n|p|exp 

4 



(m TJ f)), 02) 



K 2 n IPl^( — )• ( 93 ) 
We can insert these formulae in eqs.(18) and (83). After the change of variables 

V=1 2 H W\ ] (94) 

we obtain from the Theorem (in cq.(18) the integrals over (^-function can be 
performed directly and only the w-integral remains) 

*t(x, P) - t'- 1 / °° dr^ exp ( - X(| p | + r) ) (95) 
h{t'~ l v/|pF)*(x - nci.nr), 

where if = ^t. 
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Calculations in the light-cone coordinates for a massive particle diffusing 
without a friction (/ = in eq.(39)) are similar. We obtain from eq.(58) that 
fi = for E = O.Then, from eqs.(88)-(90) (we change variables r = p + exp(2v)) 

*t_ (P+) = P+ / °° drE [* (p+ exp(2S)) exp(2B)<5(t - A) 

^ P+ JeM^(eM2v)p + )IoC-^f 1 )5(t-^ P+ y) 

exp(-^(l + exp(2v)))dvy- 1 dy = tZ 1 /dr*(r) exp(-^)/ (f V^PT), 

(96) 

where i_ = \k 2 X-. 

We can generalize this result to a particle in a constant electric field. Then, 
from eq.(58) after a rescaling similar to eq.(87) we obtain 

^=—- ( 97 ) 

Hence, in the transition function we shall have I fl (j-^/p + r). The solution of 
eq.(73) obtained from eq.(74) is 

*tAP+) =tZ' f dr(^r*(r)eM-^)I»(±V^)- (98) 

We shall check the results (95)- (98) of the probabilistic averaging by solving 
the transport equations (73) and (81) directly. Let p = |p| then from eq.(78) p 
satisfies the stochastic equation 

3k 2 c , c 2 Bk 2 , _ i „ 
dp = —^-dt — pdt + y/cKpzdb (99) 

whereas 

dx j — — c\p\~ pjdt = —crijdt. (100) 
Eq.(lOO) has the solution 

x(i)=x-nci. (101) 

The transport equation (64) for massless particles with the Jiittner friction 
can be expressed in the form 

2 

fit* = ^-rd 2 * + ^r)0 r tf - nV x * (102) 

where a = 6k 2 c and v = [5c 2 n 2 . The transport equation (102) is related to the 
(generalized) Bessel diffusion 

K 2 C (7 K 2 C 

dtX = — d lx + — d p x - vpd P x- (103) 

The diffusion process pt solving the diffusion equation (103) (i.e., Xt — E[x(Pt)}) 
is the solution of the stochastic equation 

2 

dp= a dt- vpdt + Kdb. (104) 
2p 
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Let 

r = - A p 2 (105) 

then 

dr = ^dt - V -rdt + n^fc^db. (106) 

If v = then the Kolmogorov transition function ( the probability density to go 
from x to y in time t) for the Bessel process is (ref.[29],sec.IV.8; see also [31] [32]) 

P(t,x,y) = ±(xy) 1 -h*- 1 I i _ 1 (f ) (107) 

where a = -f- and t' = ^j^. Eq.(107) coming from ref.[29] coincides with our 
results (95)- K (9 C 6). 

We can calculate the Laplace transform of the process p+(t) applying the 
results (91), (96) and (98) of the Theorem of sec. 2 to * = cxp(— Ar) and using 
the integrals 6.614 from rcf. [33] 

£;[exp(-Ap + (i_,p + ))]=t: 1 /drexp(-Ar)exp(-^)7o(4 L N /rpT) 

- (1 + At-)- 1 cxp ( - Ap+(1 + Ai.)" 1 ) . (1 ° 8) 

The diffusion process (73) in an electric field E has the Laplace transform 

(1 + Xt-)- 1 - 2 "^' 2 exp ( - Ap+(1 + At-)" 1 ) . (109) 

The Laplace transform (108)-(109) has been calculated earlier in [29] (sec. IV. 8). 

The solution of the transport equation with friction cannot be calculated 
with Yor's formula (88). The diffusion process (see eqs.(117)-(118) below) de- 
pends on the variables A and B but in eq.(18) we still need its integral over 
time. Nevertheless, we can solve the diffusion equations (72) and (75) directly 
applying their relation to the (generalized) Bessel diffusion (104) with v > 0. 
We have calculated the transition function of the process p t with v ^ in [35] 

P t (x,y) = exp(tv(a + l))(smhvty 1 y a+1 x- a I a (isxy(smhvty 1 ) , . 
exp(-f (coth(^) + l)y 2 - f (coth(izt) - l))a; 2 ), 1 ' 

where 

v ' 2 

(a — 2 in the model of the diffusion of massless particles). The transition 
function (108) describes the imaginary time evolution in quantum mechanics 
with the potential V{x) = fx 2 + gx~ 2 (then a = \^J\ + Eg, [35]). 

The Laplace transform of the transition function of the process rt (with 
v > 0) starting from k which is the square of the generalized Bessel process can 
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be calculated using the transition function (110) (it has been has been derived 
earlier in [29]). We obtain 

£[exp(-Ar t )] = / dyP t (k, y) exp(-Ay 2 ) = ( - A K 2 ^- 1 (exp(-^) - 1) + l)~~ 
exp ( - (- AK 2 ^ 1 (cxp(-§i) - 1) + 1) Aexp(-ff)fc). 

(1U) 

For the massless particles a = 6ck 2 , v = n 2 cf3] for the process p + with the 
Jiittner friction (/' = —1) we have a = 2ck 2 and v — k 2 c/3 . The expansion of 
eq.(lll) for the diffusion of massless particles 

If °° v 

£[exp(~Ar t )] - -(c/3) 3 / dyy 2 exp(-c/3y - \y) = ^ c n {k) exp(-n-t) (112) 

n=l 

shows that dp = ^(c(3) 3 y 2 exp(—c(3y) is the invariant measure for the process 
p t and that the convergence to the equilibrium is exponential with the speed ^ 
proportional to the temperature. 

Applying the transition function (110) we can calculate the solution of the 
transport equation (81) as an expectation value over the stochastic process 

/>oo 

*t(x,|p|,n) =E[*(x-cnt,pt,n)] = / dyP t (y/\p\, ^^(x-cnt, y, n) (113) 

Jo 

From eq.(113) it follows that if "J/ as the function of x has no limit at infinity 
then the limit t — > does not exist (this is a consequence of the wave propagation). 
If the limit t — > oo in cq.(113) exists then it follows that $> t _ — » 1 and <j> t _ — » 
(up to a normalization constant). Moreover, for the momentum probability 
distribution of an "observable" <j> in a "state" <£> we obtain 



Um t ^ooJ dxdp0(p)$ t (x,p) 

= (/dp$ £ (p)0(p)) 5 L(/3c) 3 /dpdxexp(- C /3|p|)*(x,p) (U4) 

As an example of the solution (113) we could consider the plane wave with 
the initial condition ^(x, |p|,n) = cxp(i|p|nx). Then, from eq.(lll) 

*t(x,|p|,n) - (^(nx- C t)(cxp(-I« 2 c/3i)-l) + l)~ 3 

cxp ^j^(nx - ct)(cxp(-^n 2 cf3t) — 1) + 1^ exp(-i/t 2 c/3t)|p|(nx - ct)j 

(115) 

As a result we obtain again a wave (note that |p|(x — nct)n = p^x^ in the 
relativistic notation) moving in the direction n with a decreasing amplitude 
and exponentially decreasing wave vector. 
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We can solve the stochastic equations (68) and (75) and calculate expectation 
values of functions of the processes with a friction. The solution of eq.(68) reads 

|p T | = |p| exp(« 2 r + nb T ) \1 + — f3\p\ J ds cxp(n 2 s + Kb s )J (H6) 

The stochastic equation (75) for p + can also be solved with the result 

V+{t) = V+cxp{K,b T )[l + —(ip + J dscxp(Kbs)) (117) 

Unfortunately, with the Yor's formula we are unable to calculate the expectation 
value (18). However, we can obtain expectation values E[F(p T )] of any function 
of the process. This may be sufficient for a calculation of some mean values with 
the probability distribution <!> as a solution of the transport equation (11), for 
example applying eq.(18) 

J dxdp<j>{p)$(x, p) = J rfp0(p)$ £ (p)^ drE[J dx*(x,p T (p))p (T)]. 

The rhs is of the Yor's form (88) also for the processes (116)-(117) with a friction. 



8 Transport equations in a moving frame 

We consider a covariant form of the equilibrium distribution 

$ s =Po 1 exp(/(c/3%)). (118) 

The four- vector (3^ can be related to the velocity of the frame of reference in such 
a way that in the rest frame (3 = {-gf, 0, 0, 0), where T is the temperature and k 
is the Boltzmann constant (see [36] [37] [38] [39] [40] for a discussion of such equi- 
librium distributions). Then, the probability distribution da(x,p,[3) — dxdp<&E 
transforms in a covariant way under Lorentz transformations A 

da(Ax,Ap,A/3) = da(x,p, (3). 

An expectation value in a "state" $ satisfying the transport equation (11) is 

/rfprfx^(x,p)$ t (x,p) = /Ar(a:,p,^)^(x,p)S[*(x T (x,p),p T (p))]. (119) 

We applied the formula (18) in order to express the time evolution of the solution 
$ t of the transport equation by the proper time. In eq.(119) r(t) is expressed 
by the coordinate time from the equation t = (mc)" 1 pods . 
If m > then the drift (31) reads 

K j = -?^0 j f' + ^p»p j f. (120) 
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In the massless case 



ck 2 



K 3 = —^P^f- (121) 



The covariant form of the transport equation takes the form 



= £P;3§7* + + m " 2c " 2 ^)4^* (122) 



- mK c (3 + IrRVr,.. f'W. JL\D 4- gL £_JS ft,. J? 
Po ^ 

In the massless case 

(123) 

In the proper time formalism the diffusion process solving the diffusion equation 
(15) satisfies the equation 

dPj = — Pjdr —Pjf'dr + —PpMf'dT + neidb n . (124) 

where e J n (the square root of the metric in the second order differential operator 
(21), see [10]) has been defined in eq.(69). The diffusion process p r (124) is the 
same as the process p T generated by A (eq.(7)) because A = ^A^e = A. We 
have obtained this equality by direct calculations but there should be a deeper 
reason for it. 

We can derive from eq.(124) the formula for dp . Let us define 

7T = /3% (125) 

and 

X = ^x„. (126) 

Then, a direct calculation (using the Ito stochastic calculus [29]) leads to the 
proper time stochastic equations 

dn = ^. wdT+ cjf 7r 2 fdT 

-^/3"/3,frfr + ^(PoPj - PjPoHdb" 1 ; 

»-*. (128, 

In the massless case eq.(127) reads 

dir = ^firdr + ^f^fdr + nixdb. (129) 

Using the variables {X, n) we can write a proper time diffusion equation 
for an evolution of the probability distribution (13) &(X, ir) assuming that it 
depends solely on the variables (X, ir) 



20 



In the massless case this diffusion equation reads 

= irdx^ + {^ir + ^TT 2 f)d^ + ^Tr 2 d^ (131) 

The transport equation for massless particles is the same as the one for the 
Bcsscl diffusion 

d X * = yTT^* + + ^TT/'R*. (132) 

In order to preserve an interpretation of (3 as the temperature the four-vector 
(3 should be time-like. We can see from eq.(130) that the model simplifies 
substantially (the transport equation reduces to the Bessel diffusion) if we let 
— > 0. In a formal limit [3 = (/3-,0, 0, 0) (when [3 is on the light cone) and 
with 9 T $ = we obtain the model (72) corresponding to the light cone coor- 
dinates with the exception of the factor ^-d^ in cq.(130) which is replaced by 

\d+ in eq.(72). The incorrect factor in the formal limit comes from the incor- 
rect probability measure in this limit which in the light-cone coordinates should 
be dp\dpidp+p~^ exp/ and not dppy 1 exp/. The factor p^ 1 is responsible for 

the difference between the term ^f-c^ in eq.(130) and 4r<9+ in eq.(72). 

Rclativistic invariance means that a Lorentz transformation of the processes 
X and 7r can be shifted into a transformation of the frame of reference, i.e., for 
a function F 



E[F(Ax(t, (3, Ax, Ap), Ap(r, (3, Ap))] = E[F(x(t, A~ 1 f3, x, P ), P (t, A^p^p))] 

(133) 

9 Discussion 

In our earlier papers we have developed a formalism of relativistic diffusions 
evolving in the proper time. Such a formalism is useful because the explicit 
relativistic invariance is a strong guiding principle when building relativistic 
diffusion equations. In this paper we have shown that a solution of the stochastic 
equations evolving in the proper time allows a construction of the solution of 
the transport equation (in the laboratory time) by means of an integration 
over the proper time. We have shown how the method works in some models 
which can be solved exactly. These examples are not typical for the proper time 
equations. For massless particles the proper time is just an affine parameter 
on the trajectory without a physical meaning. A function depending solely 
on p + can describe a massive particle whose motion is restricted to a straight 
line (in the direction of the third axis). The results concerning the transport 
equations described by the Bessel diffusion of momenta could find applications 
in astrophysics, in plasma physics and in heavy ion collisions. The momentum 
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|p| evolves in time as the Bessel process which has some distinguished features 
among all diffusion processes [31] [32]. In particular, the (generalized) Bessel 
diffusion has the exceptional property that an exponential initial distribution 
remains exponential for any time. 

In general, the method of proper time could be applied for approximate 
calculations and computer simulations. If m > then an expansion parameter 
(mc) _1 appears in all our equations. We could use it in the Theorem (eq.(18)) 
in order to calculate the integral over r with an expanded argument of the 
8- function (p s = b s in the lowest order of the expansion in (mc)" 1 ) 

We can now calculate the expectation value (18) in a way similar to the one in 
this paper because the probability distribution of the integral dsh^. is known. 

The original relativistic diffusion in the proper time of Schay [4] and Dudley 
[5] is explicitly Lorentz invariant because its proper time evolution is generated 
by a Lorentz invariant differential operator. In applications we need diffusions 
which have a limit for a large time. The necessary drag terms which force the 
process to an equilibrium cannot be Lorentz invariant. The reason is that the 
notion of the equilibrium is itself frame dependent. We described this frame 
dependence in the last section. The stochastic process transforms in a covariant 
way with respect to the Lorentz transformation if together with the process we 
transform also the frame. 
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